Determinations of upper critical field in the c-a plane 
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Within a continuous Ginzburg-Landau model for layered superconductors, two procedures are 
proposed to determine the upper critical field parallel to the c-a crystal plane with an angle 9 tilted 
from the c-axis. For an intrinsically layered superconductor, the upper critical fields for 6 — > 90°, as 
determined by the two procedures, are consistent with each other and in reasonably good agreement 
with those determined by similar procedures suited for the parallel upper critical field {6 = 90°). 
The profile of the order parameter obtained at Bc2 is Gaussian-like, indicating the plausibility of 
the procedures proposed. 



I. INTRODUCTION 

Theories within the Ginzburg-Landau (GL) framework 
are important for studies of properties of layered super- 
conductors. Considering the layered structure of high Tc 
superconductors (HTSs), Koyama et alB first proposed 
a continuous Ginzburg-Landau (CGL) model, in which 
the GL coefficients and the superpair masses are spatially 
dependent. Recently, we have proposed a similar CGL 
model and have appliecLit to a layered superconductor 
without a magnetic fieldEl, and to layered systems with a 
magnetic field parallel to the a- (b-) crystal axiaJ. Details 
of the determinations of the parallel upper critical field 



were given in Ref. ^. In this work, our CGL formula- 
tion is applied to layered superconductors immersed in a 
magnetic field parallel to the c-a crystal plane. Two pro- 
cedures shall be proposed to determine the upper critical 
field, which is a subject of long-term interestO. 



II. MODEL 



We consider a layered superconductor comprising of 
alternating superconducting (S) and insulating (I) layers 
along the c-axis. The CGL free energy for the system iso. 
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where r — {x,y) is the planar vector and A{f,z) — 
(A^^) (r, z), y4^(r, z)) is the vector potential, ^'(f, z) is 
the superconducting order parameter and _B(r,jzlis the 
internal magnetic field. (3 is assumed constantoQ. The 
CGL condensation coefficient a{T, z) and the perpendic- 
ular and parallel afiiective masses, M(z) and m(z), are 
assumed as beforeau: 



a(T, z) = [ao + ai cos{2ttz/D)] (1 - T/Tc), 
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where ao, ai. Go, Gi, 50 and gi are the model parame- 



ters. D is the size of the unit cell equal to d//2 -f ^5 + 
di/2 = di + ds, with dj and ds denoting the thickness 
of the I and S layers, respectively. In the present simula- 
tion for layered superconductors, the intrinsically layered 
cuprate Bi2Sr2CaCu208 (Bi2212) is chosen as the mod- 
eling prototype as this superconductor has |an explicitly 
layered structure due to its high anisotropyQ. Details of 
the calculation inputs for Bi2212 can be found in Ref. ||. 

Let us now consider the case where an external mag- 
netic field B is applied in the direction tilted from the 
z-axis (c-axis) by an angle 9 in the z-x {c-a) plane. Tak- 
ing the vector potential as A = (0, B{x cos 9 — z sm9), 0), 
the linearized CGL equation from the CGL free energy 
is obtained as follows. 
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In a new coordinate system, 
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Eq. H becomes 
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Assuming ^{x',y',z') = e*^^'$(x', z'), it follows from Eq. || that 
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with 

x'q = 



hk/{2eB). For < 6* < 90°, one may choose 



For a given temperature T, the maximum magnetic 
field B which satisfies Eq. |^ gives a point on the Bc2-T 
plot. Eq. H shall be numerically solved subject to the 
following boundary conditionsQ, 



III. NUMERICAL PROCEDURES 

A. Procedure I 

Taken into account the boundary conditions, Eq. ^can 
be represented as 
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where the column vector $ = k = 1,2, {2n — 

2) X 2n represents the discrete solutions of Eq. ^ (here 
' indicates transpose, k and n are integers). The sparse 
matrix U has the following structure, 
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where the block matrices Uj i {I ~ 1,2,3,4,5) can be expressed as foUows 
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The elements ci,i,j, C2.i,j, C25,ii- are the coefficients 
of the discretized equations of Eq. g. Considering Eq. ||, 
the ranges of the index i for Ui,^, U2,i, Ua,^, XJ^^i and 
U5,j are 4 ^ 2n - 1, 3 ^ 2n - 1, 2 2n - 1, 2 ^ 2n - 2 
and 2 — > 2n — 3, respectively. The range of the index j for 
the elements in Uz,i is from 2 to 2?! — 1. We can find the 
dimension of U as follows: each matrix of the form Ua,^ 
is of order 2n; the main diagonal consists of matrices of 
the form Ua,^ and there are 2n — 2 of them (the range of 
the index i for Ua,^ is 2 ^ 2n — 1), the dimension of U 
is thus {2n - 2)2nx {2n - 2)2n. 



The boundary conditions of Eq. ^ have been incorpo- 
rated into U: the periodic condition of Eq. ^is explicitly 
expressed in all the first rows of Ua,i (see Eq. [ij); the 
derivative condition of Eq. 7b is treated by the four-point 
difference techniques in all the last rows of Ua,^; the zero 
condition of Eq. |7c|is implemented at i = 1, 2n, which is 
the reason that the range of the index i for Ua,i is from 2 
to 2n — 1. The zero solutions on and outside the bound- 
aries at i = l,2n also lead to the structures of the first 
and last two rows in U (see Eq. ^). 

Except for the first and last rows, all the other rows 
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in JJi^i have the same structure. In the second and sec- 
ond last rows of Ui^i, the periodic property of the solu- 
tions is considered. As exemplified in the second last 
rows of U;^2ri-i (l = 1,2,3,4,5), the five coefficients 

C5,2n-l,2n-l, Cio,2n-l,2n~l , Ci5.2n- 1.2rt-l , C20.2Tt-1.2ra-l 

and C25,2n-i,2n-i are there (cf. Eqs. |l^^J) because the 
discrete solutions corresponding to these coefficients are 
equivalent to the counterparts at j = 2. 

For non-trivial solutions, the determinant of U should 
be zero. 



det lUl = 0. 
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By eliminating the constant elements in the first and last 
rows (columns) in the matrices Ua^i, namely, eliminating 
the ffist and last rows (columns) in XJi^i, Eq. 15 can be 
transformed into 
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As a matter of convenience, the notations l,i axe omit- 
ted in the elements above and only the index j (j = 
2, 3, 2n — 1) is given. The quantities of v{l)i_j, 
e{l)i_j, <y(l)i_j and T{l)ij indicate the corresponding ele- 
ments from U. For example, ej in U3 ^ is e(3)ij '^i3,i,j 
in Us^j; /14 in U'^ ^ is m(1)i,4 ^ ci,i,4 in Ui,^; /14 in ^ is 
At(2)j,4 ^ C6,i,4 in U2,i; Ai4 in Ug ^ is ^(3)^^ cii,,;^ in 
Ua^i. Each Uj j is of order 2n — 2. Thus, the dimension 
of U' is (2?! - 2)(2n - 2) x (2n - 2)(2n - 2), which is less 
than that of U. 
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where e' is short for e'{l)ij and one has 
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It can be verified that the parameters of the magnetic 
field B, contained in the cis.^ j coefficients, only appear in 
the main diagonal position of U (thus U'). Since ci^^ij 
can be categorized into two different parts: with and 
without B, Eq. [l^ can be further written as 
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where I is a unitary matrix. P has the same structure 
and dimension as U' and the block matrices Pi^i in P 
have the following general structure, 
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Thus, P/,i|;=i2,4,5 = Ji = i,2,4,5 but Ps^i ^ Ugj. 

From Eq. |8| it is clear that the largest solution for S, 
namely Bc2, can be obtained from the maximum eigen- 
value of the following eigen equation. 



r 



(21) 



where % is the eigen function of P. Having obtained 
Bc2, the corresponding order parameter can be obtained 
by substituting Bc2 into Eq. |8[ 



B. Procedure II 



In the above procedure, the magnetic field square B^ 
has been treated as an eigenvalue (Eq. ^). In fact, one 
can directly discretize Eq. ^ into a matrix eigen equa- 
tion, from which the upper critical field can be directly 
deduced. 
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where the matrix Q has the same structure as U and the 
block matrices in Q have the following general structure, 
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Here, the quantities of /ij, i/j, e^- 
meaning as those in Pj.i. Unlike that in U;^i, Uj j and 
Pj^i, the range of the index j in Q; ^ is now from 1 to 
2n. The dimension of Q is the same as that of U, i.e., 
(2n — 2)2n x (2n — 2)2n. Note that the periodic prop- 
erty of the solutions has been considered implicitly in the 
first and last rows of Q; , , unlike the explicit condition in 
Ua^i (see Eq. Such a periodic property also appears 
in the second and last second rows in Q; ^ as well as in 

The present procedure allows a direct determination 
of the upper critical field from the largest eigenvalue of 
the magnetic field square of the linear GL equation, com- 
pared to the conventional method that Bc2 is obtained 
via the loweslfiiK^value of the Landau level of the linear 
GL equationtlE 



IV. RESULTS AND DISCUSSION 

In the above procedures, we have treated the magnetic 
field square as eigenvalue problems (Eq. |T] and Eq. |2^ ) 
and hence, the upper critical field in the c-a plane can be 
directly obtained from the corresponding eigen equations. 
It should be noted that these procedures are related to 
the partial differential Eq. ^, namely a 2D differential 
equation. Note that we can have similar procedures t| 
determine the upper critical field parallel to the layer: 
but these procedures (corresponding to Eqs. 5 and 6 in 
Ref. ^, respectively) are now related to an ordinary differ- 
ential equation, i.e., a ID differential equation (see Eq. 1 
in Ref. ||). The values of Bc2 of Bi2212, calculated from 
the two 2D procedures under the conditions of 6* = 89.9° 
and T = K, are presented in Table |. The correspond- 
ing ID calculations at 6* = 90° and T = K are also 
listed. It can be seen that the results obtained from the 
2D procedures {0 = 89.9°) are in good agreement with 
each other and reasonably approximate the correspond- 
ing ID calculations [6 = 90°). 

It is found that the order parameter obtained at Bc2 
from procedure I (Eq. |2l|) is also consistent with that 
from procedure II (Eq. p2[) and we shall utilize proce- 
dure II for the following calculations. In Fig. 1, we show 
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the surface and contour plots of the order parameter for 
Bi2212 with the conditions 9 = 89.9° and T = K. Note 
that the x'-axis nearly overlaps the z- (c-) axis (namely, 
x' K, z) a.i 9 = 89.9°. Hence, as depicted in Fig. 1, 
the order parameter of Bi2212 is localized in a thin slab 
(z e [—2, 2] a.u.) along the c- axis. The behavior found 
here is consistent with the fact that for Bi2212 the^sxis 
coherence length at zero temperature is very shortEj and 
less than the distance between the two effective super- 
conducting layers (~ 12.34 A, see Ref. However, for 
another intrinsically layered, superconductor YBCO with 
relatively small anisotropyllJ, we found that the localized 
domain at zero temperature is about z G [—10,10] a.u., 
which is broader than the distance between the two ef- 
fective superconducting layers (~ 8.41 A, see Ref. ||). As 
such, one may say that at zero temperature, Bi2212 ex- 
hibits a 2D feature while YBCO demonstrates a quasi-3D 
behavior. 

Fig. 2 shows the order parameter and its contour at i?c2 
for Bi2212 at a higher temperature 0.9 Tc, with the angle 
9 = 89.9° unchanged. It is obvious that the order param- 
eter becomes broader than the case at zero temperature 
(Fig. 1), as expected. The localized domain across the 
layer is now about z G [—6, 6] a.u. and it is still smaller 
than the distance between the adjacent superconduct- 
ing layers. This behavior is qualitatively consistent with 
the 2D feature of Bi2212 at 0.9 Tc calculated by the one- 
dimensional procedure. However, for YBCO, we find that 
the localized domain at 0.9 Tc is about z E [—15, 15] a.u., 
larger again than the distance between the nearest super- 
conducting layers. Consequently, one may conclude that 
Bi2212 is a 2D superconductor in a large temperature 
while YBCO is intrinsically a quasi-3D superconductor. 
Note that the order parameter presented in Figs. 1 and 
2 are found to be of a Gaussian type, which is consistent 
with the fact that the profile of the order parameter as- 
sociated with the smallest eigenvalue (Landau level) of 
the usual linear GL equation at Bc2- 

Finally, it is worth mentioning that Eq. 2l| of pro- 
cedure I requires less memory than Eq. |22 of proce- 
dure II in obtaining Bc2 since the dimension of P is 
less than that of Q. The number of the stored ele- 
ments in the sparse matrix P is less than that in Q by 
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(2n-2)2nx(2n-2)2n-(2n-2)(2n-2)x(2n-2)(2n-2) = 
16(2n — l)(n — 1)^. When n is large, the reduced amount 
(~ 3271"^) is significant. However, in order to obtain both 
Bc2 and the associated order parameter, procedure I re- 
quires more CPU time than procedure II since in using 
procedure I, two equations (Eqs. ^ and have to be 
solved while only one equation (Eq. p2[ ) is involved in 
procedure II. 



V. CONCLUSION 

In this paper, we applied our continuous Ginzburg- 
Landau model to layered superconductors in the presence 
of a magnetic field parallel to the c-a plane. By treating 
the magnetic field square as eigenvalues, two procedures 
were proposed to determine the upper critical field. Near 
the direction parallel to the layers, the critical fields of 
Bi2212 calculated from the two procedures are consis- 
tent with each other and are good approximations to the 
corresponding parallel upper critical fields determined by 
similar procedures. The behaviors of the order parameter 
obtained at i?c2 a-re reasonable. Our procedures can be a 
useful starting point for investigating the properties asso- 
ciated with the upper critical field and the corresponding 
order parameter. 
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TABLE L Values of 5^2 (Tesla) at 89.9° at K for Bi2212, 
determined by the two 2D procedures (Eq. ^ and Eq. p2[ ). 
The results from the 2D procedures are consistent with each 
other and reasonably approximate the parallel upper critical 
fields determined by the corresponding ID procedures. 



procedure 
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Bc2 


I (2D, Eq. §) 


30 


3342.3065 


40 


3337.3166 


50 


3336.7392 


I (ID) 


800 


3338.6584 




30 


3338.7414 


II (2D, Eq. H) 


40 


3336.9549 




50 
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H (ID) 


800 


3338.6801 



FIG. 1. Surface and contour plots of the order parameter 
of Bi2212 at zero temperature. The plots are depicted against 
the z' — x' plane {x' is from negative to positive and z' is from 
5000 to 12000). The order parameter is localized in a narrow 
domain (about 6. a.u.) along the c-axis unit cell. 



FIG. 2. Surface and contour plots of the order parameter 
of Bi2212 at 0.9 Tc- The order parameter spreads out and 
becomes broader than that at zero temperature (see Fig. 1). 
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